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O+HCl  Chemistry  Models  for  Hypervelocity  Collisions 

in  DSMC 
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and 

I.  J.  Wysong* 

Edwards  Air  Force  Base,  AFRL-PRSA,  Edwards  AFB,  CA  93524 


In  earlier  work,  atmospheric- jet  interaction  flows  were  simulated  using  direct  simulation 
Monte  Carlo  (DSMC)  calculations  for  altitudes  of  80,  120  and  160  km.  At  high  altitudes, 
the  0(3P)+HC1(* 1S+)  — >  OH(2n)+Cl(2P)  reaction  was  found  to  contribute  the  most  to  OH 
production.  However,  the  total  collision  energy  (TCE)  reaction  probability  was  often  found 
to  be  unphysically  greater  than  unity.  In  this  work,  we  examine  in  detail  different  reaction 
models  such  as  the  TCE  model  using  the  rate  constants  of  Mahmud  et  al  and  Xie  et  al.  In 
addition  to  the  TCE  model,  quasi-classical  trajectory  (QCT)  calculations  were  performed 
to  obtain  the  reaction  probability  for  the  O+HCl  reaction  with  the  new  benchmark  triplet 
A”  and  A’  surfaces.  Both  reaction  and  total  collision  cross  sections  were  calculated  by  the 
QCT  method.  The  dynamic  molecule  collision  model  was  used  to  calculate  the  viscosity 
cross  section  from  which  the  VHS-equivalent  collision  cross  sections  were  derived.  The 
tabulated  QCT  reaction  probabilities  were  then  used  in  the  DSMC  calculations  to  model 
OH  production  at  120  km.  It  is  found  that  the  QCT-based  collision  cross  sections  are 
greater  than  the  Bird  VHS  cross  sections,  and  the  QCT  reaction  probabilities  for  O+HCl 
are  lower  than  the  TCE  probabilities  using  the  rate  of  Mahmud  et  al.  Using  the  QCT 
reaction  probability,  the  maximum  reaction  probability  for  the  O+HCl  reaction  was  found 
to  be  lower  than  0.4  in  the  atmospheric  -  jet  interaction  flows  for  an  altitude  of  120  km 
and  a  freestream  velocity  of  5  km/s. 


Nomenclature 


X  Deflection  angle 

rj  Viscosity  coefficient 

e'ini  Dimensionless  post-collisional  internal  energy 

Hint  Dimensionless  internal  energy 

etr  Dimensionless  translational  energy 

g  Dimensionless  relative  velocity 
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g'  Dimensionless  post-collisional  relative  velocity 
fj,  Reduced  mass 

v  Coefficient  of  viscosity 

ui  Viscosity  index  (w  =  u  —  0.5) 

(f>  Initial  direction  of  rotational  vector 

ip  Initial  orientation  of  molecules 

o>  Reaction  cross  section 

ox  Total  collision  cross  section 

Viscosity  cross  section 
<jei  Elastic  collision  cross  section 

< Ttr  Inelastic  collision  cross  section 

< 7j-hs  VHS  collision  cross  section 

Cv  Effective  number  of  degrees  of  freedom  for  vibration 

b  Impact  parameter 

d  Diameter 

Ec  Collision  energy 

Ea  Activation  energy 

Eint  Internal  energy 

Etr  Translational  energy 

g  Relative  velocity 

J  Rotational  quantum  number 

k  Boltzmann  constant 

Kf  Reaction  rate  constant 

to  Mass  of  a  molecule 

Pr  Reaction  probability 

Qint  Partition  function  for  internal  modes 

Qrot  Rotational  partition  function 

QVib  Vibrational  partition  function 

v  Vibrational  quantum  number 


I.  Introduction 

The  modeling  and  simulation  of  chemically  reacting  flows  caused  by  the  interaction  of  a  jet  positioned 
on  the  side  of  a  hypersonic  rocket  is  a  challenging.  The  forward  and  aft  reaction  control  systems  (RCS) 
engines  provide  the  thrust  for  attitude  maneuvers  to  control  the  vehicle  trajectory.  The  investigation  of  the 
interaction  of  RCS  jets  with  the  rarefied  atmosphere  is  needed  at  high  altitudes  since  the  relative  velocities  of 
the  freestream  atomic  oxygen-plume  chemical  species  are  sufficiently  high  to  produce  chemical  species  that 
can  radiate  in  various  portions  of  the  optical  spectrum.  We  consider  a  generic  RCS-vehicle  geometry  with 
freestream  conditions  corresponding  to  those  of  high  altitude,  hypersonic  flight  for  a  lateral  side  jet  thrusting 
perpendicular  to  the  rocket  velocity  vector  into  a  near-vacuum  environment.  An  important  feature  of  the 
jet  interaction  are  the  hypervelocity  chemical  reactions  between  the  jet  and  atmospheric  species  that  occur 
for  reactant  conditions  very  far  from  the  conditions  for  which  the  Arrhenius  kinetic  rates  were  obtained. 

In  earlier  work,1  we  demonstrated  the  utility  of  direct  simulation  Monte  Carlo  (DSMC)5  simulations 
at  80,  120,  and  160  km  altitudes.  The  entire  set  of  the  chemical  reactions  between  thruster  side  jet  and 
plume-atmospheric  species  includes  reactions  between  oxygen  and  nitrogen  species  as  well  as  ones  that 
either  produce  or  consume  OH  (see  Table  1).  The  complete  list  of  chemical  reactions  modelled  in  the  DSMC 
simulations  may  be  found  in  Ref.  1.  The  first  three  reactions  in  Table  1  involve  the  dissociation  of  water  by 
freestream  constituents,  N2,  02,  and  O.  The  relative  importance  of  these  three  reactions  will  change  with 
altitude.  The  last  four  reactions  are  exchange  reactions  between  the  freestream  and  plume  constituents.  The 
following  two  exchange  reactions,  0+H2  and  O+HCl,  are  potentially  important  because  both  H2  and  HC1 
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are  present  in  relatively  large  mole  fractions  of  divert  solid  propellant  motors.  At  altitudes  of  120  km  and 
higher,  O  exchange  with  HC1  was  found  to  contribute  more  than  70  %  of  the  OH  produced.1  Also,  it  was  not 
possible  to  establish  a  consistent  trend  in  terms  of  the  relative  contributions  of  the  OH  producing  chemical 
reactions  as  a  function  of  altitude. 

A  detailed  study  of  the  DSMC  flowfield  calculation  at  120  km  altitude  for  a  freestream  velocity  of  5  km/s 
illustrated  an  important  problem  in  the  reaction  model  for  O+HCl.  A  schematic  of  the  flow  geometry  for 
the  RCS  jet  plume-atmospheric  interaction  is  given  in  Fig.  1.  Figure  2  shows  the  spatial  distribution  of  the 
total  collision  energy  (TCE)6  reaction  probability  for  O+HCl  — >  OH+C1  in  the  X-Y  plane.  The  interaction 
region  ahead  of  the  jet  is  seen  for  an  altitude  of  120  km  and  a  freestream  velocity  of  5  km/s.  It  is  observed 
that  the  TCE  reaction  probabilities  for  this  reaction  are  greater  than  one  in  the  interaction  region.  Figure 
3  shows  the  distribution  of  the  TCE  reaction  probability  for  all  of  the  O+HCl  — >  OH+C1  reactions,  with 
the  result  that  more  than  10  percent  of  the  reactions  have  a  probability  greater  than  one.  Since  the  reaction 
probability  cannot  be  greater  than  one,  the  specific  implementation  used  in  Ref.  1  was  to  artificially  limit  the 
reaction  in  all  cases.  This  essentially  causes  a  much  smaller  number  of  reactions  to  occur  in  the  simulation 
compared  to  that  governed  by  the  Arrhenius  equation. 

For  hypervelocity  collisions,  such  as  occur  in  the  atmospheric  -  jet  interaction  flows,  the  extension  of 
the  reaction  rate  for  the  O+HCl  reaction  to  temperatures  higher  than  3,000  K  and  the  subsequent  use  of 
the  rate  in  the  TCE  model  are  problematic.  In  order  to  reevaluate  the  adequacy  of  the  present  model  for 
the  O+HCl  reaction,  it  was  decided  to  improve  upon  the  TCE  model  by  use  of  a  more  exact,  fundamental 
approach  that  uses  the  molecular  dynamics/quasi-classical  trajectory  (MD/QCT)  method  to  calculate  the 
reaction  collision  cross  sections.  In  addition,  the  Dynamic  Molecular  Collision  (DMC)  model  was  used  to 
calculate  both  total  collision  cross  sections  and  viscosity  collision  cross  sections,  and  the  MD/QCT  viscosity 
collision  cross  sections  were  compared  with  those  of  the  variable  hard  sphere  (VHS)  model. 

Table  1.  Freest  ream-plume  species  reactions  for  OH  production  reactions  used  in  the  TCE  model 


Reaction 

A,  m3/s 

n 

Ea,  x  10-19  J 

h2o+n2->oh+h+n2 

5.81  X  10-15 

0.00 

7.314 

h2o+o2->oh+h+o2 

1.13  X  10-7 

-1.31 

8.197 

h2o+o->oh+h+o 

1.13  X  10-7 

-1.31 

8.197 

h2o+o->oh+oh 

1.13  X  10-16 

0.00 

1.275 

h+o2->oh+o 

1.66  X  10-16 

0.00 

1.061 

o+h2->oh+h 

3.12  x  10-16 

0.00 

0.952 

OH+C1— s-O+HCl 

3.10  X  10-27 

2.91 

0.070 

O+HCl— s-OH+Cl(Mahmud)2 

5.60  X  10-27 

2.87 

0.244 

O+HCl— s-OH+Cl(Xie)3 

1.70  X  10-22 

1.485 

0.408 

The  MD/QCT  modeling  builds  on  the  work  of  Ozawa  et  al'  wherein  the  TCE  chemical  reaction  model6 
was  compared  with  MD/QCT  calculations  for  reactions  related  to  the  production  of  OH  in  a  bow-shock. 
The  MD/QCT  method  enables  the  physical  investigation  into  reaction  cross  sections,  viscosity  cross  sections 
as  well  as  total  collision  cross  sections.  It  was  found  that  at  higher  altitudes,  the  discrepancy  between  the 
TCE  and  MD  model  was  larger  due  to  the  greater  thermal  non-equilibrium.  The  MD/QCT  model  is  less 
affected  by  the  translational  energy,  compared  to  the  TCE  model,  but  it  is  more  sensitive  to  the  reactant 
internal  energy. 

The  outline  of  the  remainder  of  this  paper  is  as  follows.  The  MD/QCT  modeling  of  the  reaction  and 
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total  cross  sections  is  presented  in  Sec.  II  and  compared  with  the  results  of  the  TCE  model  using  newer 
Arrhenius  parameters  as  well  as  more  accurate  O+HCl  viscosity  values  at  high  temperatures.  In  Sec.  Ill,  we 
discuss  the  specific  DSMC  numerical  parameters  used  for  the  single  freestream  condition  that  we  examine 
here,  120  km,  5  km/s.  Finally,  in  Sec.  IV,  we  present  detailed  comparisons  of  the  MD/QCT  results  with  the 
TCE  model  and  compare  the  changes  observed  in  the  DSMC  simulation  for  the  different  chemical  reaction 
parameters  for  the  O+HCl  reaction.  Although  there  are  multiple  chemical  reactions  that  contribute  to  the 
production  of  OH,  we  limit  on  discussion  to  the  O+HCl  reaction  since  this  is  the  most  energetically  favorable 
reaction  path. 

II.  MD/QCT  Modeling  of  the  Reaction  Probability  and  Total  Cross  Section 

A.  MD/QCT  Reaction  Probability 

The  widely  used  method  for  calculating  the  reaction  probabilities  in  DSMC  is  the  TCE  model.6,6  In  this 
model,  the  reaction  probability  has  a  special  form  that  allows  one  to  match  experimental  reaction  rates 
Kf(T)  in  modified  Arrhenius  form, 

Kf  =  ATn  ejqy(—Ea/kT),  (1) 

where  A  is  the  preexponential  factor,  n  is  the  temperature  dependence,  and  Ea  is  the  activation  energy. 
The  Arrhenius  data  presented  in  the  Table  1  were  used  to  calculate  the  reaction  cross  sections  based  on  the 
TCE  model.  In  our  previous  work,1  the  rate  of  Mahmud  et  al ,2  A  =  5.6  x  10-27  m3/molec/s,  n  =  2.87,  and 
Ea  =  2.44  x  10_2°  J,  was  used  for  O+HCl— »OH+Cl,  but  in  this  work,  a  newer  Arrhenius  expression  based 
on  the  calculations  of  Xie  et  al 3  was  used. 

A  parallel  MD  code  was  developed  to  calculate  the  reaction  and  viscosity  cross  sections  using  the  QCT- 
internal  energy  quantum  mechanical  threshold  (QCT-IEQMT)  method.3  For  the  O+HCl— »OH+Cl,  the  ab 
initio  3  A”  potential  energy  surface  of  Ramachandran  and  Peterson  (RP) 1  was  utilized.  This  surface  was  also 
utilized  by  Xie  et  al,  who  computed  both  QM  and  QCT  reaction  cross  sections  for  a  range  of  energies.1  The 
present  results  agree  closely  with  those  of  Xie’s  for  the  same  energy  conditions.  However,  it  was  necessary 
for  us  to  perform  our  own  MD/QCT  calculations  since  a  large  set  of  probabilities  are  required  for  the  DSMC 
flowfield  simulation,  and  we  require  cross  sections  for  total  collision  energy  values  of  Ec  up  to  8  x  10-19  J. 
In  addition,  DSMC  uses  continuum  vibrational  and  rotational  energies,  requiring  a  table  of  probabilities  for 
a  range  of  Eint ,  while  Xie  et  al  computed  reaction  cross  sections  for  v=0,l,  and  2  states  only. 

The  MD/QCT  reaction  probability  is  the  ratio  of  the  reaction  cross  section  obtained  by  the  QCT- 
IEQMT  method  to  the  total  collision  cross  section.  An  efficient  Monte  Carlo  numerical  procedure  is  utilized 
to  evaluate  the  multidimensional  integrals  associated  with  averaging  of  the  scattering  geometric  properties 
before  a  collision.  The  Monte  Carlo  method  converges  at  a  rate  that  is  independent  of  the  dimensionality  of 
the  integral.  In  order  to  utilize  the  Monte  Carlo  method,  the  specific  initial  states  of  each  of  the  trajectories 
that  correspond  to  the  desired  reaction  cross  section,  such  as  <r(T),  a(v,j,Etr)  or  a{Eint,Etr)  are  needed. 
Microcanonical  sampling  enables  one  to  efficiently  sample  the  initial  states  of  the  target  molecule  of  equal 
energy  to  obtain  the  trajectory  initial  conditions  necessary  to  perform  the  MD/QCT  calculations.10  For  the 
DSMC  simulations,  the  reaction  probability  of  a  chemical  reaction  as  a  function  of  translational  energy  Etr 
and  molecular  internal  energy  Eint,  is  used.  The  reaction  cross  section  is  obtained  by  evaluating  whether  a 
O+HCl  collision  ends  in  a  chemical  reaction  by 


(J  r 


=  Tib 


2 

max 


Nr 

At 


(2) 


where  Nr  is  the  number  of  trajectories  that  result  in  a  reaction  and  Nt  is  the  number  of  total  trajectories. 
In  previous  work, 1  the  reaction  probability  was  assumed  to  be 


Pr  = 


(J f 

&VHS 


(3) 
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where  the  collision  cross  section  was  found  to  be  very  close  to  the  VHS  cross  section  using  Bird  values.5  For 
the  O+HCl  case,  however,  it  will  be  shown  that  the  MD/QCT  total  cross  section  is  different  than  that  given 
by  the  VHS  model  and  a  different  relationship  will  be  used  to  obtain  the  MD/QCT  reaction  probability. 


B.  Total  Collision  Cross  Section 

As  mentioned  above,  in  our  previous  work,  the  variable  hard  sphere  (VHS)  model  with  Bird  values  was  used 
for  the  total  collision  cross  section.  The  VHS  cross  section  is  defined  by  a^HS  =  tt d2  with6 

d  =  dref[2kTref  / (fxg2)w /T(2  -  a;)]1/2,  (4) 

where  to  =  0.375,  dref  =  4.38  A,  and  Tref  =  273  K  for  the  O+HCl  reaction  in  Ref.  5.  The  reference 
diameters  are  calculated  from  the  viscosity  data, 


j  _  ,  30(mkTref  /n)l/2  ^1/2 

Tef  ~  M(5  -  2w)(7-  20j)fj,ref  ’ 

One  could  also  attempt  to  calculate  the  total  collision  cross  section  for  the  RP  potential  surface,  but  the 
unbounded  nature  of  a  classical  total  cross  section  for  a  realistic  potential  comes  into  play.  To  illustrate,  one 
may  use  the  MD/QCT  method  to  compute, 


<?t  =  nb2maxr^-  (6) 

where  Nc  is  the  number  of  trajectories  that  result  in  a  collision  and  Nt  is  the  number  of  total  trajectories. 
However,  in  order  to  decide  whether  a  collision  has  occurred,  the  cut-off  deflection  angle  x  is  needed,  and 
<jt  was  found  to  be  strongly  dependent  on  the  cut-off  angle  x  for  the  O+HCl  system. 

Figures  4  and  5  show  the  distribution  of  the  deflection  angle  x  for  O+HCl  collision  at  Eint  =  0.5  x  10-19 
and  2.0  x  10-19  J,  respectively,  for  three  different  relative  collision  velocities  obtained  from  our  MD/QCT 
calculations.  These  HC1  internal  energies  were  chosen  because  these  numbers  are  in  the  energy  range  of 
interest  in  the  flowfield  freestream  conditions  at  120  km  and  velocity  of  5  km/s.  As  may  be  seen  in  the 
figures,  the  translational  and  internal  energies  increase,  the  peak  deflection  angle  x  shifts  to  small  angles. 
Note  that  the  deflection  angle  distribution  shown  in  Fig.  4  is  based  on  a  smaller  range  of  impact  parameters 
than  that  for  Fig.  5.  Our  QCT  procedure  only  considered  trajectories  out  to  the  impact  parameter  value 
where  the  reaction  probability  became  negligible.  Thus,  the  higher  case  considers  trajectory  samples 
from  higher  impact  parameter  values.  This  fact  leads  to  the  visibly  lower  occurrence  of  high  deflection  angles 
for  the  Eint  =  2  x  10-19  J  case  in  Fig.  5  for  v  =  3  km/s.  The  MD/QCT  collision  cross  section  results 
change  with  different  selected  cut-off  angles.  (The  question  of  whether  different  values  of  internal  energy 
should  be  assigned  different  DSMC  collision  cross  sections  has  been  discussed,11  but  will  not  be  considered 
here.)  Thus,  for  higher  internal  and  translational  energies,  the  cut-off  angle  is  more  important,  and  it  is 
difficult  to  define.  Therefore,  the  DMC  method  is  utilized  to  calculate  the  viscosity  cross  section. 

For  implementation  in  DSMC,  it  is  more  consistent  to  obtain  a  transport-based  collision  cross  section 
(usually  viscosity)  from  which  to  convert  a  reaction  cross  section  to  a  reaction  probability,  since  the  DSMC 
method  already  uses  the  VHS  model  viscosity-based  cross  section  as  the  basis  for  computing  the  number  of 
collisions  per  cell  per  time-step  and  as  the  basis  for  collision  dynamics  (hard-sphere  isotropic  scattering).  In 
this  spirit,  the  work  of  Tokumasu  and  Matsumoto12  demonstrated  the  Dynamic  Molecular  Collision  (DMC) 
model  to  calculate  an  accurate  viscosity  cross  section  for  the  given  potential  using  a  Monte  Carlo  method  for 
integration.  While  the  integration  over  the  impact  parameter  becomes  infinite  for  the  total  collision  cross 
section,  those  for  classical  momentum  and  energy  transfer  cross  sections  are  finite.  Therefore,  the  viscosity 
cross  section  was  calculated  first,  and  converted  to  the  equivalent  VHS  collision  cross  section  for  each  collision 
velocity.12 
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The  viscosity  cross  section  is  calculated  by  the  Monte  Carlo  evaluation  of  an  integral  given  by  Ref.  13 


9  ■  2  |  1  (  A  ~  \2 

-j  sin  x+-(Aei„t) 


-(Aei„t)2sin2  x)<1t. 


(7) 


where  g  is  defined  as  \/ [x/kTg ,  e,„t  is  nondimensionalized  by  kT ,  and  Ae,„t  is  the  change  of  the  dimensionless 
HC1  internal  energy  before  and  after  a  collision.  The  quantities  x>  9  and  Ae,„t  are  obtained  from  the 
similar  MD/QCT  calculations  for  the  reaction  probability  except  the  sampling  of  initial  trajectory  conditions 
corresponds  to  the  geometric  pre-collisional  conditions  specified  by  Ae,„t.  The  MD  viscosity  cross  section 
converges  when  the  maximum  impact  parameter  is  selected  so  that  the  effect  of  potential  is  neglected. 
The  integral  /()dr  specifies 


/<  >*-r  u o  €&*«**« 

where  bmax  is  the  cut-off  impact  parameter.  In  addition,  the  viscosity  coefficient  rj  is  obtained  by 
^  g  i  r  °°  r°°  r°°  tp.  2 


Tj  5\Z27T flkT  Qint 


1  C°°  r°°  C° 

—  dv  dJ 

int  Jo  Jo  Jo 


1 7-1  Etr  (n  j  ,  1  \  (  ^9  Evj  \ 


(8) 

(9) 


The  equivalent  VHS  collision  cross  section  is  ut  =  7TG?f ,  and  the  diameter  for  the  total  collision  cross  section 
can  be  obtained  by 


dt  = 


6  a. 


H,MD 


7T3 


(10) 


III.  DSMC  Models,  Numerical  Parameters,  and  Freestream  Conditions 

Figure  1  shows  the  flow  geometry  used  in  the  modelling  of  the  interaction  of  the  atmosphere  and  the  RCS 
jet.  A  small  rocket  is  modeled  as  a  blunted  cone  cylinder  and  a  thruster  positioned  on  the  cylinder  right 
after  the  cone-cylinder  junction.  The  radius  of  the  cylinder  is  0.2  m,  and  the  length  from  the  head  of  the 
cone  to  the  nozzle  exit  is  2  m.  In  this  work,  the  angle  of  attack  is  zero.  The  freestream  parameters  at  120  km 
altitude  are  listed  in  Table  2.  A  starting  surface  was  obtained  from  the  axisymmetric  DSMC  plume  core-flow 
DSMC  simulations.  The  plume  core-flow  DSMC  simulations  were  performed  with  the  nonuniform  nozzle 
exit  condition.1’15  The  density  isolines  of  about  6  x  1021  molecules/m3  were  taken  for  the  starting  surface 
of  a  60-lbf  (270-N)  thruster.  The  starting  surface  is  an  oval  shape  with  approximate  x  and  y  dimensions  of 
0.3  and  0.5  m. 


Table  2.  Freestream  parameters 


Parameter 

120  km 

Temperature,  K 

354 

Number  density,  molec/m3 

4.73  x  1017 

O2  mole  fraction,  % 

9 

N2  mole  fraction,  % 

73 

O  mole  fraction,  % 

18 

The  three-dimensional  DSMC  calculations  were  implemented  in  the  SMILE  computational  tool.  The 
gas  is  considered  a  14-species  reacting  mixture.  Approximately  18  million  molecules  were  simulated  in 
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the  computational  domain.  Separate  grids  were  used  for  collisions  and  macroparameters  adaptive  to  flow 
gradients.  The  total  number  of  collision  cells  was  approximately  5  million.  The  total  number  of  time  steps 
was  about  100,000  with  a  time  step  of  2.0  x  10-7  s.  Macroparameter  sampling  was  started  after  20,000  time 
steps,  a  time  sufficient  to  reach  the  steady  state.  These  parameters  were  chosen  to  minimize  the  statistical 
dependence  between  simulated  particles,  remove  the  grid  dependence  of  the  results,  and  furnish  sufficient 
spatial  resolution  of  the  boundary  layer  along  the  rocket. 

The  majorant  frequency  scheme16  is  employed  for  modeling  the  molecular  collision  frequency  to  enhance 
the  statistical  representation.  For  trace  species,  a  weighting  factor  of  approximately  0.01  was  used.  The 
VHS  model  was  used  for  modeling  non-reactive  interactions  between  particles.  The  Borgnakke-Larsen 1 
(BL)  model  (also  see  Ref.  18)  with  temperature-dependent  rotational  and  vibrational  relaxation  numbers 
was  chosen  for  modeling  rotation  -  translation  and  vibration  -  translation  energy  transfer.  The  effective 
degree  of  freedom  for  vibration  was  modeled  to  be  temperature-dependent.  The  TCE  model  was  used  to 
calculate  the  chemical  reaction  probabilities.  For  the  O+HCl— »OH+Cl  reaction,  both  the  TCE  and  MD 
models  were  used.  The  entire  set  of  the  chemical  reactions  can  be  found  in  earlier  work,1  and  the  OH 
production  reactions  used  in  this  work  are  presented  in  Table  1.  The  gas-surface  interaction  was  modeled 
using  a  diffuse  model  with  total  energy  and  momentum  accommodation  with  a  rocket  wall  temperature  of 
300  K. 


IV.  Results  and  Discussion 

A.  Rate  Constant  for  O+HCl  — >  OH+C1 

Although  the  DSMC  method  requires  the  probability  of  a  reaction,  instead  of  the  rate  constant,  the  latter 
quantity  provides  a  consistency  check  for  the  MD/QCT  calculations  presented  in  this  work.  Figure  6  presents 
the  reaction  rate  constant  calculated  using  a  parallel  MD  code  between  1,000  K  and  3,000  K.  The  reaction  rate 
may  be  calculated  from  the  reaction  cross  section  by  averaging  the  target  molecule  (HC1)  over  a  Maxwellian 
equilibrium  distribution, 


kf 


f(T)QvibQrot 


1/2 


rOO  rbmax 

VV(2J+1)  /  Pr(v,J,Etr,b ) 

„  T  J0  Jo 


V  J 


x  exp(— £+  / kT^nbdbEtrdEtr i 


(11) 


and  for  O+HCl, 


f(T)  = 


3 

5  +  3  exp(— 228/T)  +  exp(-326/T) ' 


(12) 


The  factor  f(T)  is  the  probability  that  the  target  molecular  system  is  initially  on  one  of  the  three  electronic 
surfaces  that  allow  a  reaction  to  occur.  This  temperature  dependent  expression  accounts  for  the  spin- 
orbit  splitting  of  the  overall  triplet  reagents.  The  typical  number  of  trajectories  per  temperature  value  was 
approximately  20,000.  Figure  6  shows  that  the  rate  of  Mahmud  et  al 2  used  in  our  previous  DSMC  calculations 
is  higher  than  the  QCT(RP)  rate  obtained  with  the  RP  surface  in  this  work  and  the  improved  canonical 
variational  theory  (ICVT)  rate  constant  of  Xie  et  al.  The  QCT  (RP)  rate  is  slightly  lower  than  the  ICVT 
rate  because  only  the  3  A"  surface  was  used  for  QCT  calculations,  and  the  3 A'  surface  contribution  increases 
for  higher  temperature.  In  addition,  the  tunneling  effect,  not  included  in  the  MD/QCT  calculations,  is  more 
important  for  lower  temperature. 


7  of  18 


American  Institute  of  Aeronautics  and  Astronautics  Paper  2006-1193 
Approved  for  public  release;  distribution  is  unlimited 


B.  Reaction  Probability  for  O+HCl  — >  OH-j-Cl 

Figure  7  shows  a  comparison  of  reaction  probabilities  between  the  MD  and  TCE9  models  as  a  function  of 
reactant  internal  energy  at  3  and  5  km/s  relative  velocities.  The  MD  reaction  probability  is  calculated  here 
as  the  reaction  cross  section  divided  by  the  VHS  cross  section  using  Bird’s  values.  -  The  rate  of  Mahmud  et 
al 1  used  for  our  previous  simulations  and  the  rate  of  Xie  et  al 3  are  also  shown  in  the  figure.  The  effective 
number  of  degrees  of  freedom  for  vibration  Cv,hci  =  1  is  used,  and  the  VHS  model  with  the  parameters  of 
Ref.  5  of  w  =  0.375,  dref  =  4.38  A  and  Tref  =  273  K  was  used  for  all  the  three  cases  in  this  subsection. 
In  the  next  subsection,  the  fidelity  of  the  VHS  values  used  will  be  discussed.  It  can  be  seen  that  the  TCE 
probabilities  computed  from  the  rate  of  Mahmud  et  al  are  much  higher  than  the  MD/QCT  or  TCE  model 
using  the  Arrhenius  parameters  derived  from  the  rate  coefficients  computed  by  Xie  et  al  for  both  relative 
velocities.  This  result  is  consistent  with  our  previous  calculations  in  which  the  probabilities  for  the  O+HCl 
reaction  were  found  to  be  greater  than  one.  For  the  3  km/s  case,  the  MD/QCT  model  predicts  slightly  higher 
reaction  probabilities  than  the  TCE  with  the  rate  calculations  of  Xie  et  al.  Nonetheless,  at  5  km/s,  the 
discrepancy  becomes  more  significant,  and  the  MD/QCT  model  predicts  higher  probabilities  than  the  TCE 
model.  The  reason  for  this  is  as  follows.  For  the  TCE  model  the  reaction  probability  increases  uniformly 
as  the  translational  energy  increases.  However,  the  MD  reaction  probability  does  not  necessarily  increase 
uniformly  as  the  translational  energy  increases. 

C.  Collision  Cross  Section  for  O+HCl 

In  the  previous  subsection,  the  VHS  collision  cross  section  with  ui  =  0.375,  dref  =  4.38  A  and  Tref  =  273 
K  was  used.  In  this  subsection,  the  preferred  values  for  VHS  for  O+HCl  collisions  is  investigated.  In  our 
earlier  reported  results,1’15  the  viscosity  index  lu  of  0.25  for  O  and  u>  of  0.5  for  HC1  were  used,5’19  i.e.,  the 
HC1  molecule  was  treated  in  the  Maxwell  model.  However,  although  the  viscosity  index  oj  of  HC1  is  0.5 
between  20  and  99  °C,  for  some  gases,  w  decreases  as  temperature  increases.19  Therefore,  the  VHS  total 
cross  sections  obtained  from  a  Maxwell  molecule  may  not  be  accurate  for  higher  temperatures. 

Let  us  consider  the  selection  of  the  two  parameters  of  coefficient  of  viscosity  u  and  the  reference  diameter 
dref  for  the  VHS  model.  The  parameters  used  in  the  previous  DSMC  calculations  were  mostly  obtained  from 
viscosity  data  in  the  low  temperature  range.19  Figure  8  shows  the  variation  of  viscosity  coefficient  of  HC1 
from  100  to  5,000  K.  In  the  figure,  the  data  of  Svehla20  and  the  data  (Bird)  in  Ref.  19  are  shown.  Svehla 
estimated  the  viscosities  using  the  Lennard-Jones  (12-6)  potentials.  From  the  data  (Bird),  the  coefficient 
of  viscosity  of  HC1  vhci  is  1.0,  and  with  vhci  of  1.0,  the  coefficient  of  viscosity  is  higher  than  the  data  of 
Svehla.  The  coefficient  of  viscosity  of  HC1  vhci  of  0.65  agrees  well  with  the  data  of  Svehla  between  1,000  and 
5,000  K.  The  VHS  parameters  for  HC1  are  listed  in  Table  3.  It  is  seen  that  the  value  of  w  is  quite  different 
for  low  and  high  temperatures. 


Table  3.  VHS  parameters  for  HC1 


Data 

dref 5  ^ 

Tref ,  K 

Bird 

0.50 

5.76 

273 

Svehla 

0.15 

3.15 

5,000 

Using  the  RP  3  A"  surface  in  the  MD/QCT  calculations,  the  quantities  such  as  deflection  angle  and  the 
change  in  HC1  internal  energy  are  obtained  and  used  in  Eq.  (7)  to  calculate  the  viscosity  cross  section. 
The  change  in  HC1  internal  energy,  ej„t,  is  obtained  by  analyzing  the  momentum  of  the  post-collisional 
HC1  species  in  each  trajectory.  To  account  for  the  neglect  of  the  RP  3 A"  potential,  the  maximum  impact 
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parameter  was  selected  such  that  the  rate  of  translational  energy  transfer  was  less  than  10  %,  and  the  MD 
viscosity  cross  section  o^md  converged.  The  DMC  method  is  used  in  this  work  to  obtain  both  the  O+HCl 
viscosity  coefficients  as  well  as  the  VHS-equivalent  collision  cross  section. 

The  O+HCl  viscosity  coefficient  for  high  temperatures  was  calculated  by  the  DMC  method.  Figure  9 
shows  a  comparison  of  the  viscosity  coefficients  obtained  from  the  MD/QCT  results  and  the  O  and  HC1  data 
in  Ref.  19  (Bird)  and  data  of  Svehla."-  To  compare  with  the  MD/QCT  result  we  take  the  average  of  the 
separate  O  and  HC1  viscosity  coefficients.  The  MD/QCT  results  for  O+HCl  are  believed  to  be  the  most 
accurate  and  fall  above  or  below  the  data  of  Refs.  19  (Bird)  and  20,  depending  on  the  temperature.  The 
agreement  between  the  MD/QCT  results  and  the  viscosities  derived  from  the  low-temperature  data  of  Ref. 
19  is  good.  At  higher  temperatures,  the  difference  in  the  MD/QCT  results  and  the  viscosities  derived  from 
the  coefficients  of  Svehla  is  attributed  to  the  simpler  potential  model  used  in  the  latter  calculations. 

Figure  10  presents  a  comparison  of  the  collision  cross  sections  between  the  DMC(MD/QCT)  and  VHS 
models  at  Eint  =  0.5  x  10-19  and  2.0  x  10-19  J.  For  the  two  viscosity  parameters  tested  in  the  VHS  model, 
Fig.  10  shows  that  the  collision  cross  sections  are  greater  with  high  temperature  data  of  Svehla.2  The 
VHS  cross  section  with  ljhci  =  0.15  and  duel, ref  =  3.15  at  5,000  K  (Svehla)  is  significantly  larger  than 
the  VHS  cross  section  with  ujhci  =  0.5  and  duel, ref  —  5.76  at  273  K  (Bird)  by  more  than  50  %.  The 
MD/QCT  results  are  between  the  VHS  cross  sections  obtained  from  the  data  of  Ref.  5  (Bird)  and  those 
of  Svehla.20  Therefore,  if  the  MD/QCT  VHS-equivalent  collision  cross  section  in  Eq.  (10)  is  used,  lower 
reaction  probabilities  are  predicted  than  those  that  would  be  obtained  from  the  VHS  cross  sections  using  the 
data  in  Ref.  5  (Bird).  Also,  it  was  found  that  the  MD/QCT  VHS-equivalent  cross  section  does  not  change 
significantly  with  changes  of  the  HC1  internal  energy.  Since  the  MD/QCT  VHS-equivalent  cross  section  is 
not  dependent  on  the  HC1  internal  energy,  the  cross  sections  could  be  fit  to  the  simpler  VHS  model.  The 
parameters,  ui  —  0.39  and  dre/,ioooif  =  3.9  A  were  found  to  give  a  good  fit  of  the  MD/QCT  cross  sections 
to  the  VHS  form.  The  total  cross  sections  derived  from  these  parameters  will  be  designated  as  ■ 

The  inelastic  collision  cross  section  is  defined  for  collisions  that  do  not  react  but  produce  a  change  in  Etr 
as, 

atr  =  J  (Aetrfdr,  (13) 

and  is  related  to  the  total  collision  cross  section  as 

Or  =  CHr  +  crei  +  0>.  (14) 

The  quantity  Aetr  is  the  change  of  the  dimensionless  translational  energy  before  and  after  a  collision  and  was 
calculated  from  the  MD/QCT  pre  and  post  trajectory  HC1  conjugate  momenta  for  non-reactive  trajectories 
as, 

o'2  52 

Aetr  =  —  +  ecM  — y  (15) 

where  ecM  is  the  dimensionless  center-of-mass  translational  energy.  Figure  11  shows  the  inelastic  collision 
cross  section  as  a  function  of  the  relative  velocity.  As  shown  in  the  figure,  the  inelastic  cross  section  decreases 
as  the  translational  energy  increases.  It  is  also  found  that  the  inelastic  cross  section  is  a  function  of  both 
the  translational  and  internal  energy.  However,  because  the  change  of  atr  is  insignificant  compared  to  the 
VHS-equivalent  collision  cross  section,  a^^p  does  not  show  a  strong  dependence  on  internal  energy. 

Figures  12  and  13  show  a  comparison  of  the  reaction  probability  for  O+HCl  -+  OH+C1  calculated  by 
using  the  reaction  cross  section  divided  by  either  the  VHS  model  (Bird)  or  at  3  km/s  and  5  km/s, 

respectively.  The  VHS  model  with  w  =  0.39  and  dre/,ioooif  =  3.9  A  was  used  to  represent  a^^fD  because 
the  MD/QCT  collision  cross  section  results  agree  well  with  the  VHS  model  with  those  parameters  as  shown 
in  Fig.  10.  Both  the  MD/QCT  and  TCE  reaction  probabilities  using  a^^fD  are  lower  than  those  with 
the  VHS  cross  section  (Bird,  uihci  =  0.5)  because  a^^fD  are  greater  than  the  VHS  cross  sections  (Bird, 
(jJhci  —  0.5).  However,  the  change  in  the  TCE  model  is  smaller  than  that  in  the  MD/QCT  model. 
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Figure  14  presents  the  distributions  of  the  Aej„i  at  5  km/s  for  O+HCl  collision  for  two  initial  internal 
energies.  A imt  was  calculated  by  A e*nt  =  e'int  —  imt ■  Figure  14  shows  that  as  the  initial  HC1  internal 
energy  increases,  the  A e,„t  distribution  spreads.  In  other  words,  the  inelastic  cross  section  increases,  as  the 
initial  internal  energy  increases,  as  is  generally  expected  based  on  simple  models  (see  e.g.,  Ref.  21).  At 
Eint  =  0.5  x  10-19  J,  energy  is  mostly  transferred  from  the  translational  to  HC1  internal  modes.  Figure  15 
shows  the  distributions  of  the  A e»„t  at  Eint  =  2.0  x  10-19  J  for  different  collisional  translational  energies. 
At  3  km/s,  the  energy  is  transferred  from  the  HC1  internal  to  translational  mode.  However,  as  the  initial 
translational  energy  increases,  the  energy  shifts  from  the  translational  to  internal  modes.  In  general  the 
results  show  a  tendency  toward  energy  equilibration  between  modes.  Where  Etr  >  Eint.  the  Eint  gains  and 
where  Eint  >  Etr.  the  Etr  gains. 

D.  DSMC  Calculations 

First,  the  TCE  model  with  the  rate  constant  of  Xie  et  al  for  the  O+HCl  -+  OH+C1  reaction  was  utilized  for 
the  DSMC  calculation  at  120  km  altitude  for  freestream  velocity  of  5  km/s.  For  the  TCE  model,  the  VHS 
total  collision  cross  section  for  O+HCl  with  ui  =  0.375,  dref  =  4.38  A  and  Tref  =  273  K  was  used.  Figure  16 
shows  the  N2  number  density  (molecule/m3 )  contour  for  this  case.  Note  that  the  change  of  chemical  reaction 
models  does  not  affect  the  overall  flow  field.  The  maximum  N2  number  density  is  about  1.2  x  1018  m-3 
between  1  and  2  m  ahead  of  the  rocket  cone.  Figure  17  shows  the  translational  temperature  contour  with 
the  rate  constant  of  Xie  et  al  for  the  O+HCl  -+  OH+C1  reaction  at  120  km  altitude  for  freestream  velocity  of 
5  km/s.  Similar  to  the  total  number  density,  the  overall  translational  temperature  profile  does  change  with 
the  different  reaction  probabilities  for  the  O+HCl  reaction.  While  the  high  number  density  shock  region  is 
located  between  x  =  -2  and  -1  m,  the  high  temperature  shock  region  is  shifted  more  forward  than  the  high 
number  density  shock  region.  Figure  18  shows  the  Mach  number  contour  at  120  km  altitude  for  freestream 
velocity  of  5  km/s.  The  high  temperature  shock  region  and  low  Mach  number  region  are  nearly  the  same. 

Secondly,  the  DSMC  calculations  were  implemented  with  three  chemical  reaction  models  at  120  km 
altitude  for  freestream  velocity  of  5  km/s.  For  the  O+HCl  — >  OH+C1  reaction,  the  TCE  model  with  the 
previous  rate  constant2  (1)  and  the  rate  constant  of  Xie  et  al 3  (2)  and  the  MD/QCT  reaction  probability 
(3)  were  used.  The  MD/QCT  reaction  probability  for  O+HCl  — >  OH+C1  used  in  this  subsection  is  the 
ratio  of  the  MD/QCT  reaction  cross  section  to  Vtmdi  and  MD/QCT  reaction  probabilities  were  tabulated 
to  be  used  in  DSMC.  Figure  19  shows  the  OH  number  density  contours  with  the  TCE  model  of  the  rate 
(2)  for  O+HCl  at  120  km  altitude  for  freestream  velocity  of  5  km/s.  The  maximum  OH  number  density 
was  about  2  x  1016  m-3  with  this  chemical  reaction  rate.  In  Figs.  20  and  21,  the  comparison  between 
the  three  chemical  reaction  models  is  shown.  Figure  20  presents  the  OH  number  density  distribution  along 
the  y  =  8  m  line.  Figure  21  shows  the  OH  number  density  distribution  along  the  x  =  —1.5  m  line.  As 
expected,  the  TCE  model  predicts  higher  OH  production  than  the  MD/QCT  model.  With  the  rate  (1),  the 
OH  was  produced  much  more  than  the  other  two  cases  (2)  and  (3)  because  the  probability  for  O+HCl  was 
significantly  high  as  shown  in  Fig.  2  and  7.  With  the  rate  (2),  the  maximum  OH  number  density  is  reduced 
by  more  than  50  %,  and  with  the  MD/QCT  probability,  it  is  reduced  by  a  factor  of  4.  In  the  shock  region, 
the  reaction  probabilities  changed  dramatically  between  the  rate  (1)  and  (2)  for  the  TCE  model.  With  the 
rate  (1),  the  reaction  probabilities  result  in  greater  than  one.  (In  the  DSMC  calculation,  they  are  set  to 
one.)  Due  to  the  jet-atmosphere  interaction,  high  translational  energies  result  in  high  reaction  probabilities 
for  the  TCE  model.  About  10  percent  of  the  O+HCl  collisions  results  in  the  probability  greater  than  one 
with  the  rate  (1)  using  the  TCE  model  as  shown  in  Fig.  3.  In  contrast,  with  the  rate  of  Xie  et  al  (2), 
most  of  the  O+HCl  collisions  are  predicted  to  have  the  reaction  probability  smaller  than  0.4  (see  Fig.  22). 
The  reaction  probability  with  the  MD/QCT  is  further  lower  than  the  other  two,  and  more  than  80  %  of 
the  reaction  had  the  probability  of  lower  than  0.1  as  shown  in  Fig.  23.  Using  the  same  VHS  (Bird  values) 
collision  cross  section  between  the  TCE  and  MD/QCT,  the  TCE  with  the  rate  (2)  predicts  higher  reaction 
probability  than  the  MD/QCT  as  shown  in  Fig.  7.  In  contrast,  if  both  the  reaction  and  total  collision  cross 
sections  are  obtained  by  the  MD/QCT  method,  the  MD/QCT  reaction  probability  becomes  lower  than  the 
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TCE  with  the  rate  (2).  At  5km/s,  the  probabilities  between  the  MD/QCT  and  the  TCE  with  the  rate  (2) 
are  very  close,  but  at  3km/s  or  8km/s,  the  TCE  probabilities  with  the  rate  (2)  predicts  higher  than  those 
with  the  MD/QCT. 


V.  Conclusion 

The  chemical  reacting  flows  with  the  interaction  of  a  jet  have  been  simulated  using  the  DSMC  method. 
The  objective  of  this  work  was  to  improve  the  fidelity  of  the  chemical  reaction  models  in  DSMC,  and 
analyze  the  most  important  chemical  reaction,  0(3P)+HC1(1E+)  -»  0H(2II)+C1(2P)  in  the  atmosphere 
-  jet  interaction  flows.  The  MD/QCT  calculations  were  performed  for  the  O+HCl  reaction  using  the  ab 
initio  3 A"  potential  energy  surface  of  Ramachandran  and  Peterson  (RP)1  for  both  reaction  and  collision 
cross  sections.  The  DMC  model  was  used  to  calculate  the  viscosity  cross  sections,  and  the  VHS-equivalent 
collision  cross  sections  were  obtained  with  the  assumption  of  isotropic  scattering.  Finally,  the  MD/QCT 
reaction  probabilities  were  compared  with  those  obtained  by  the  TCE  model  with  two  rates. 

These  results  show  that  using  the  same  VHS  collision  cross  section,  the  MD/QCT  reaction  probability 
was  predicted  to  be  lower  than  the  TCE  probability  with  the  rate  of  Mahmud  et  al2  used  in  Ref.  1  (1)  and 
higher  than  the  TCE  probability  with  the  rate  of  Xie  et  al 3  (2).  The  TCE  model  with  the  rate  (1)  predicted 
the  reaction  probability  higher  than  one  in  the  energy  range  of  interest.  In  general,  the  TCE  model  provides 
reasonable  results  if  it  uses  reasonable  Arrhenius  pararmeter  values  and  if  there  is  no  strong  favoring  in  the 
reaction  among  different  forms  of  reactant  energy. 

The  VHS  cross  section  is  very  sensitive  to  the  viscosity  parameters,  and  there  is  something  of  uncertainty 
in  it.  Therefore,  to  be  more  accurate,  the  MD/QCT  total  collision  cross  sections  using  the  DMC  model  were 
calculated  for  each  condition.  The  energy  transfer  between  the  translational  and  internal  modes  was  also 
investigated.  While  the  translational  energy  transferred  to  the  internal  energy  at  the  initial  HC1  internal 
energy  of  5  x  10-20  J,  the  average  of  energy  shift  was  almost  zero  at  the  initial  HC1  internal  energy  of 
2  x  10-19  J.  However,  the  magnitudes  of  the  inelastic  cross  sections  are  small  compared  to  the  total  cross 
section.  The  dependence  on  the  initial  internal  energy  is  so  small  that  the  MD/QCT  VHS-equivalent  collision 
cross  sections  were  converted  to  fit  the  VHS  with  uj  =  0.39  and  dref  =  3.9  at  1,000  K.  These  values  were 
used  for  the  MD/QCT  reaction  probabilities,  and  these  probabilities  were  tabulated  to  be  used  in  DSMC. 

The  DSMC  simulations  were  performed  for  the  atmosphere  -  jet  interaction  flows  for  an  altitude  of  120 
km  and  a  freestream  velocity  5  km/s.  The  chemical  reaction  models  did  not  affect  the  overall  flowfield,  but 
affected  the  properties  of  product  species,  OH.  Although  the  TCE  with  the  rate  (1)  resulted  in  the  reaction 
probability  higher  than  one,  the  TCE  with  the  rate  (2)  and  MD/QCT  models  showed  the  O+HCl  reaction 
with  the  probability  less  than  0.4.  Also,  the  OH  production  was  slightly  lower  with  the  MD/QCT  than  the 
TCE  with  the  rate  (2). 
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Figure  1.  Schematic  of  the  flow. 


Figure  2.  Reaction  probability  distribution  in  xy- 
plane  for  O+HCl  — >•  OH+C1  using  the  TCE  model 
(rate  constant  of  Mahmud  et  al 2)  for  an  altitude  of 
120  km  and  a  freestream  velocity  5  km/s. 
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Figure  3.  Reaction  probability  distribution  for 
O+HCl  — >■  OH+C1  using  the  TCE  model  (rate  con¬ 
stant  of  Mahmud  et  al)  for  an  altitude  of  120  km  and 
a  freestream  velocity  5  km/s. 


Figure  4.  Distribution  of  the  deflection  angle  x  f°r 
O+HCl  collision  at  Eint  —  0.5  X  10-19  J:  Dependence 
on  the  translational  energy. 
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Figure  5.  Distribution  of  the  deflection  angle  x  f°r 
O+HCl  collision  at  Eint  =  2.0  X  10-19  J:  Dependence 
on  the  translational  energy. 


Figure  7.  Comparison  between  MD  and  TCE  mod¬ 
els:  Reaction  probabilities  for  O+HCl  — >•  OH+C1  as  a 
function  of  the  reactant  internal  energy. 


Figure  6.  Reaction  rate  constant  for  O+HCl  — >•  Figure  8.  HC1  viscosity  coefficient  for  high  tempera- 
OH+C1  as  a  function  of  temperature.  The  QCT(RP)  ture. 
is  calculated  in  this  work. 
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Figure  9.  O+HCl  viscosity  for  high  temperatures  cal¬ 
culated  by  the  MD/QCT  method  (Eq.  9). 


Figure  10.  Comparison  of  total  collision  cross  sections 
between  MD/QCT  and  VHS  models  for  O+HCl  colli¬ 
sion. 


Figure  11.  Comparison  of  inelastic  collision  cross  sec¬ 
tions  for  O+HCl  collision  at  Eint  =  0.5  X  10“ 19  and 

Eint  =  2.0  X  10-19  J. 


Figure  12.  Comparison  of  the  reaction  probability  cal¬ 
culated  by  using  the  reaction  cross  section  divided  by 
either  the  VHS  model  (Bird)  or  the  MD/QCT  total 
collision  cross  section  (<j^(MD)):  Reaction  probabili¬ 
ties  as  a  function  of  the  reactant  internal  energy  at  3 
km/s. 
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Figure  13.  Comparison  of  the  reaction  probability  cal¬ 
culated  by  using  the  reaction  cross  section  divided  by 
either  the  VHS  model  (Bird)  or  the  MD/QCT  total 
collision  cross  section  (<jt(MD)):  Reaction  probabili¬ 
ties  as  a  function  of  the  reactant  internal  energy  at  5 
km/s. 
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Figure  15.  Distribution  of  the  A e%nt  for  Eint  =  2.0  X 
10“ 19  J  at  3  km/s  (top),  5  km/s  (middle),  and  8  km/s 
(bottom)  for  O+HCl  collision. 


Figure  14.  Distribution  of  the  A£int  for  E{nt  —  0.5  X 
10“ 19  J  (top)  and  E{nt  =  2.0  X  10“ 19  J  (bottom)  at  5 
km/s  for  O+HCl  collision. 
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Figure  16.  N2  number  density  (molecule/m3)  contours  Figure  18.  Mach  number  contours  with  the  rate  of  Xie 
with  the  rate  of  Xie  et  al  (TCE)  at  120  km  altitude  for  et  al  (TCE)  at  120  km  altitude  for  freestream  velocity 
freestream  velocity  of  5  km/s.  Area  shown  is  7.4  X  7.5  of  5  km/s.  Area  shown  is  7.4  X  7.5  m. 
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Figure  17.  Translational  temperature  contours  with 
the  rate  of  Xie  et  al  (TCE)  at  120  km  altitude  for 
freestream  velocity  of  5  km/s.  Area  shown  is  7.4  x  7.5 
m. 


Figure  19.  OH  number  density  (molecule/m3)  con¬ 
tours  with  the  rate  of  Xie  et  al  (TCE)  at  120  km  alti¬ 
tude  for  freestream  velocity  of  5  km/s.  Area  shown  is 

7.4  x  7.5  m. 
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Figure  20.  Comparison  of  OH  number  density 
(molecule/m3)  at  y  =  8  m  for  three  chemical  reaction 
models  at  120  km  altitude  for  freestream  velocity  of  5 
km/s. 


Figure  21.  Comparison  of  OH  number  density 
(molecule/m3)  at  x  =  —1.5  m  for  three  chemical  reac¬ 
tion  models  at  120  km  altitude  for  freestream  velocity 
of  5  km/s. 
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Figure  22.  Reaction  probability  distribution  for 
O+HCl  — }  OH+C1  using  the  TCE  model  (the  rate  of 
Xie  et  at)  for  an  altitude  of  120  km  and  a  freestream 
velocity  5  km/s. 


Figure  23.  Reaction  probability  distribution  for 
O+HCl  — >■  OH+C1  using  the  MD  model  for  an  alti¬ 
tude  of  120  km  and  a  freestream  velocity  5  km/s. 
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